Higher-nodal collective modes in a resonantly interacting Fermi gas 
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We report on experimental investigations of longitudinal collective oscillations in a highly elon- 
gated, harmonically trapped two-component Fermi gas with resonantly tuned s-wave interactions 
('unitary Fermi gas'). We focus on higher-nodal axial modes, which in contrast to the elementary 
modes have received little attention so far. We show how these modes can be efficiently excited 
using a resonant local excitation scheme and sensitively analyzed by a Fourier transformation of 
the detected time evolution of the axial density profile. We study the temperature dependence of 
the mode frequencies across the superfluid phase transition. The behavior is qualitatively different 
from the elementary modes, where the mode frequencies are independent of the temperature as 
long as the gas stays in the hydrodynamic regime. Our results are compared to theoretical predic- 
tions based on Landau's two-fluid theory and available experimental knowledge of the equation of 
state. The comparison shows excellent agreement and thus both represents a sensitive test for the 
validity of the theoretical approach and provides an independent test of the equation of state. The 
present results obtained on modes of first-sound character represent benchmarks for the observation 
of second-sound propagation and corresponding oscillation modes. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.Lm 



I. INTRODUCTION 

In ultracold quantum gases, measurements on collec- 
tive oscillations are well established as powerful tools to 
study the many-body properties of the system [1, 2]. Ex- 
periments on collective modes reveal the dynamics in 
the different regimes of superfluid, collisionally hydro- 
dynamic, and collisionless behavior. The eigenfrequen- 
cies can be determined very accurately, which allows to 
extract valuable information on the equation of state 
(EOS), with great sensitivity to subtle interaction effects 
in the strongly interacting regime. 

In ultracold Fermi gases [3-5], collective modes have 
been widely applied to study the crossover from Bosc- 
Einstcin condensation (BEC) to a Bardeen-Cooper- 
Shrieffer (BCS) type superfluid. A situation of particular 
interest is the two-component Fermi gas with resonant 
interactions, with an s-wave scattering length tuned to 
infinity by means of a Feshbach resonance [6] . This spe- 
cial case, which lies right in the center of the BEC-BCS 
crossover, has attracted a great deal of interest, mainly 
attributed to its universal properties. The resonantly in- 
teracting Fermi gas is characterized by strong interaction 
effects in the EOS [7-10] and reveals a unique universal 
thermodynamic behavior [11]. 

So far, experiments on collective modes in harmoni- 
cally trapped Fermi gases have been restricted to a few 
elementary modes. The most simple modes, sloshing 
modes, do not provide any information on the properties 
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of the quantum gas, and their main application is thus to 
accurately determine the trap frequencies. Surface modes 
are insensitive to the EOS, but they allow to clearly dis- 
tinguish between hydrodynamic and collisionless behav- 
ior [12, 13] and they have been used to detect the angular 
momentum in a rotating Fermi gas [14, 15]. Elementary 
compression modes of axial [16, 17] and radial [16, 18- 
22] character have been very widely studied in the field. 
Such modes do not only probe the particular collision 
regime, but they also give access to the compressibility 
of the gas. However, for a unitary Fermi gas, the eigen- 
frequencies of the simple compression modes do not show 
any variation across the superfluid phase transition as the 
temperature is varied [20, 22]. This can be understood as 
a consequence of the fact that superfluid and collisional 
hydrodynamics both lead to the same frequencies. A 
rigorous proof for this temperature-independence can be 
given in terms of an exact scaling solution of the hydro- 
dynamic equations of motion [23] . The situation is strik- 
ingly different for higher-nodal modes. Here the frequen- 
cies vary across the superfluid phase transition, when the 
dynamical regime changes from superfluid to collisional 
hydrodynamics [24, 25]. Such higher- nodal modes there- 
fore represent an interesting addition to the experimental 
tool-box for probing strongly interacting Fermi gases. 

We have recently carried out a series of experiments 
on higher-nodal axial modes in the geometry of a highly 
elongated trapping potential. First results on the tem- 
perature dependence have already been presented in 
Ref. [24], and the general theoretical framework is de- 
scribed in Ref. [25] . In this Article, we briefly summarize 
the main theoretical predictions (Sec. II), we describe the 
experimental procedures in more detail (Sec. Ill), and we 
present the whole set of our experimental results obtained 
for two different higher- nodal modes (Sec. IV). While we 



here restrict our attention to modes of first-sound char- 
acter, we note that the results are important as bench- 
marks for the observation of second-sound propagation 
[26] and in view of future experiments on second-sound 
modes (Sec. V). 



II. THEORETICAL PREDICTIONS 

Higher-nodal collective modes in Fermi gases have 
been theoretically studied based on Landau's two-fluid 
equations for an isotropic harmonic trapping geometry 
[27, 28]. For real experiments, however, the situation 
of highly elongated harmonic traps is more relevant. In 
this geometry, the description can be reduced to a set of 
effectively ID hydrodynamic equations, which only de- 
pend on the axial coordinate z. This simplification leads 
to a powerful approach to describe sound propagation 
and collective modes in experimentally realistic situa- 
tion. The basic approach was introduced in Ref. [29] 
for a cylindrical trap geometry with tight radial confine- 
ment. Reference [25] presents a generalization to the sit- 
uation of additional weak axial confinement, which read- 
ily describes the commonly used geometry of a highly 
elongated trap containing a 'cigar-shaped' atomic sam- 
ple. Here, we summarize the main elements of this the- 
oretical approach and the corresponding predictions for 
higher-nodal modes of first-sound character. 

The two basic assumptions underlying the ID hydro- 
dynamic approach are thermal equilibrium in the radial 
direction and a flow field that is independent of the ra- 
dial position. This corresponds to conditions of sufficient 
heat conductivity and sufficient shear viscosity, which are 
readily satisfied for resonantly interacting Fermi gases 
under common trapping conditions. Applying the local 
density approximation, one can describe the thermody- 
namics and the flow properties of the trapped sample 
using effective ID quantities, which are derived by in- 
tegrating over the transverse degrees of freedom, such 
that a thermodynamic quantity q yields a ID counter- 
part qi = 2ir J Q qrdr. 

For a first-sound collective mode with frequency w, the 
local flow speed can be expressed as v(z, t) = v z (z)e~ ZUJt , 
where the z-dependent amplitude v z (z) represents the 
spatial oscillation of the flow velocity. The hydrodynamic 
equation that describes v z (z) takes the form [25] 
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Here, ui z represents the trap frequency along the axial 
direction, m is the atomic mass, P\ is the 'ID pressure' 
(having units of force), and n\ is the linear number den- 
sity. The equation is valid for small-amplitude oscilla- 
tions, which can be treated as perturbations by lineariz- 
ing Landau's equations. 

At zero temperature and in the classical limit of high 
temperature the hydrodynamic equation (1) admits ana- 
lytic solutions of polynomial form v z = akZ k +cik-2Z k ~ 2 + 



..., with integer values of k. At T = 0, where Pi (ni)/ni = 
(2/7) [/i — (I/2)mw 2 z 2 ], with /j, being the chemical po- 
tential at the center of the trap, the frequency of the fc-th 
mode is given by 

W 2 = i(fc + I)(fc + 5K 2 . (2) 

In the classical limit, where P\jn\ = fc^T, one finds the 
different fc-dependence 

^ 2 = i(7fc + 5K 2 . (3) 

We point out that Eqs. (2) and (3) give the same values 
for k = (center of mass oscillation) and k = 1 (lowest 
axial breathing mode). In fact, one can prove that the 
frequencies of these two lowest modes are temperature 
independent for the resonantly interacting gas (unitary 
Fermi gas) , corresponding to an exact scaling solution of 
the two-fluid hydrodynamic equation [23]. On contrary, 
the frequencies of the k > 2 modes vary with tempera- 
ture. 

We now focus on the modes with k = 2 and k — 3, 
which are experimentally most relevant. Using a varia- 
tional approach [30], one can obtain their eigenfrequen- 
cies at finite temperatures as 
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In these equations t 2 = M Q Mi/M% and t 3 = M 2 M 6 /Mf , 
where we have introduced the dimcnsionless moments 
rPvo 

dx(Pno-xf +1}/2 fn(x). (6) 



Here, the phase-space density f n {x) is a universal func- 
tion [II] defined by f n (x) = nA^, where n is the 3D 
number density and At = h/(2nmkBT) 1 / 2 is the ther- 
mal dcBroglic wavelength. The dimcnsionless parameter 
x = /3/i, with f3 = 1/ksT and /i being the chemical poten- 
tial, is related uniquely to T/T F . The universal function 
f n (x) can be determined from the recent EOS measure- 
ments [7-10]. In this work, we make use of the latest 
results from [10]. 

One can also show that the velocity fields for the k = 2 
and k — 3 modes take the form 
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Here, the parameter xq = /3/^o is the value of x at the 
center of the trap. Finally, using the equation of con- 
tinuity under the ID formulation d t ni + d z (niv z ) = 0, 
one can calculate the shape of the density oscillations for 
each mode. 



III. EXPERIMENTAL PROCEDURES 



density profile n\{z) is given by [9, 33] 



A. Sample preparation 

The starting point of our experiment is an ultracold, 
resonantly interacting Fermi gas in an elongated optical 
dipole trap. This gas is prepared by evaporating a bal- 
anced mixture of fermionic 6 Li atoms in their two low- 
est spin states at a magnetic field of 834 G, very close 
to the center of the well-known broad Feshbach reso- 
nance [6, 31]. The atomic cloud contains typically N/2 = 
1.5 x 10 5 atoms per spin state. For the lowest tempera- 
tures, the waist of the trapping beam (wavelength 1075 
nm) is 31 fim, the trap depth is about 2 fiK, and the ax- 
ial and radial trap frequencies arc cu z = 2ir x 22.52(2) Hz 
and uj r = 2tt x 473(2) Hz, respectively. For experiments 
at higher temperatures, the beam waist is increased to 
38 /zm, and deeper traps are used (up to 16 /xK depth) 
with trap frequencies of up to lo z — 2ir x 23.31(3) Hz and 
Lu r = 2ir x 1226(6) Hz. The corresponding Fermi temper- 
atures Tp = h^Nujfuiz) 1 / 3 /ks vary between about 0.8 
and 1.5 fiK. 

To achieve the lowest possible temperatures, we per- 
form deep evaporative cooling up to the point where the 
trapping potential cuts slightly into the Fermi sea, indi- 
cated by the onset of spilling losses in the last stage of 
the evaporation. After that, the gas is adiabatically re- 
comprcssed by increasing the trapping beam's power to 
the extent where the trap depth becomes at least twice 
more than the Fermi energy koTp. This recompression 
step is essential to ensure negligible anharmonicities in 
the radial confinement. The essentially perfect harmonic 
confinement along the axial direction is ensured by the 
magnetic trapping that results from the curvature of the 
magnetic field used for Feshbach tuning [32]. 

We vary the temperature T of the gas by controlled 
heating, always starting from a deeply cooled cloud 
(T/T F ra 0.1). In the low-temperature range (T < 
0.2 Tp), we simply introduce a variable hold time of up 
to 4 s in which residual trap heating slowly increases the 
cloud's temperature. For the higher-temperature range 
(0.2 Tp < T < 0.5 Tp), we heat the sample using para- 
metric heating, modulating the trap power at about 2w r , 
and introducing a sufficient hold time to reach thermal 
equilibrium between the different degrees of freedom. We 
note that we use deeper traps for samples with higher 
T /Tp, because plain evaporation puts a limit on the max- 
imum attainable temperature of the gas. 



B. Thermometry 

We determine the temperature T of the gas by analyz- 
ing its density distribution in the trap, based on knowl- 
edge of the EOS. Under the local density approximation, 
one can readily show for a harmonic trap that the ID 
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Here, the universal function f p (x) is related to the uni- 
versal function f n (x) introduced in Sec. II by f p (x) = 
J_ fn(y)dy, and is therefore also known from a given 
EOS. The parameter xq is the value of x at the center of 
the trap. 

In the ideal case where one is able to obtain an accu- 
rate in-situ measurement of n\{z) by imaging the trapped 
cloud, it is straightforward to retrieve the parameters T 
and xq by fitting ni(z) using Eq. (9). However, in real- 
ity we have to deal with imperfections of our absorption 
imaging scheme. To extract the temperature in an ac- 
curate way, we have adopted the methods described in 
detail in the Appendix. 



Exciting and observing higher-nodal collective 
modes 



We apply a resonant excitation scheme to create a col- 
lective oscillation. Figure 1(a) illustrates the basic geom- 
etry of our scheme in which a repulsive 532-nm green laser 
beam perpendicularly intercepts the trapping beam. To 
excite a mode, we position the green beam near the antin- 
ode of the mode and modulate its power at the expected 
frequency of the mode. The amplitude, duration and 
shape of the modulation are carefully adjusted in order 
not to overdrive the excitation while maintaining suffi- 
cient signal-to-noise ratio. We adopt an excitation pulse 
that contains 8 cycles of sinusoidal modulation with a 
half-cycle sine envelope [34], sec illustration in Fig. 1(b). 
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FIG. 1. (Color online) Experimental scheme to excite higher- 
nodal first-sound longitudinal modes. In (a), we illustrate 
the basic geometry of exciting the optically trapped cloud 
with a weak, power-modulated repulsive laser beam, which 
perpendicularly intersects the trapping beam. In (b), we show 
the power modulation of the repulsive beam for the excitation. 



We set the maximum barrier height of the green beam to 
about 0.01 kgTp. Depending on the order of the mode 
to be excited, this is realized with beam waists ranging 
from 30 ^m to 70 fim and values of the maximum power 
Pq between 400 /iW and 3 mW. The 8-cycle pulse is cho- 
sen such that the resulting total excitation duration is not 
too long as compared to the damping time of the highest- 
nodal mode that we can observe. The smooth half-cycle 
sine envelope reduces the Fourier width and avoids side 
lobes in the spectrum, thereby suppressing the excitation 
of unwanted modes. For an efficient excitation of a given 
mode, we find that the width of the green beam should 
well match the local mode profile at the selected antin- 
ode, while the excitation frequency should be within 1% 
of the actual mode frequency. 

Once a collective mode is excited, we record the axial 
density profiles n\ (z, t) of the gas for a variable time delay 
t after the excitation pulse, where ni(z, t) is the number 
density integrated over the transverse degrees of freedom. 
These profiles are obtained with a probe beam that per- 
pendicularly intercepts the trapping beam, and are taken 
600 /j,s after suddenly releasing the atoms from the optical 
trap [35]. To enhance the visibility, we subtract a back- 
ground profile n\{z) obtained from averaging the profiles 
over all measured delay times. This gives a differential 
density variation function 8n\(z,t) — m(z,t) — fii(z), 
which is finally normalized to the maximum linear den- 
sity ni(0) at the trap center. In the top panel of Fig. 2, 
we show examples of this signal for the k = 1, 2, and 
3 modes for the coldest samples with T/Tp « 0.1. One 
can see that the adjacent antinodes always oscillate in op- 
posite directions, similar to standing waves on a guitar 
string. 



it is possible by simply analyzing the peaks in the corre- 
sponding Fourier spectrum, we adopt the following al- 
gorithm. We project 5ni(z,t) onto the mode profile 
8n\{z) to obtain a mode amplitude function A(i) = 
/_ Sni(z,t)5fii(z) dz. Then we fit a simple damped 
harmonic oscillation to A(t) to obtain the frequency and 
the damping time of the mode. This projection proce- 
dure is analogous to the projection of a superposition 
wavefunction onto one of the orthogonal cigenstates of a 
quantum system. It greatly enhances the signal-to-noise 
ratio and results in very low statistical uncertainties for 
the mode frequencies, with relative uncertainties as low 
as in the range of a few permillc. 

In Fig. 3, we show examples of A(t) for the k = 1, 2, 
and 3 modes obtained from samples with T/Tp « 0.1. 
For the k = 1 compression mode (upper panel), the ob- 
served behavior does not show any damping. Even data 
taken after a much longer delay time of 3 s (not shown) 
do not reveal any significant damping. In contrast, the 
higher-nodal modes show clear damping. At the lowest 
temperatures, the 1/e damping time for the k = 2 mode 
(middle panel) is about 1.5 s, and that for the k = 3 mode 
(lower panel) is about 0.3 s. 

We finally note that we have not succeeded in observ- 
ing modes with k > 4, despite of considerable efforts. 
We believe that this is due to a fast increase of damping 
with the mode order, which is clearly indicated by our 
data for the k = 1,2,3 modes. Large damping affects 
both our resonant excitation scheme and the detection of 
the mode by means of a Fourier transform, which may ex- 
plain a huge difference between the observed mode with 
k = 3 and the unobserved mode with k = A. 



D. Analyzing the eigenmodes: Extracting mode 
profiles, frequencies, and damping rates 

The first step to analyze the observed time-dependent 
profiles 6ni(z,t) is a Fourier transform, yielding a repre- 
sentation of our data in frequency space, 5n\(z, w) [36]. 
For this purpose we employ a fast Fourier transform 
(FFT) algorithm. Corresponding results are shown in 
the middle panel of Fig. 2, as calculated for each time- 
dependent oscillation profile in the top panel. The dis- 
crete nature in frequency space becomes evident, with 
very little noise in the background. 

It is straightforward to extract the mode profiles 
5n\{z) from the FFT results by setting Sh'l(z) = 
Sfii(z,u)k), where Wfc is the eigenfrequency of the k th 
mode. The corresponding mode profiles for the k = 1, 2, 
and 3 modes are shown in the bottom panel of Fig. 2. Ex- 
perimentally, we make use of the mode profiles and the 
frequencies obtained in this way to optimize the beam 
waist and the modulation frequency in our excitation 
scheme. We proceed iteratively, which eventually allows 
for an optimum excitation of a single mode. 

To extract the mode frequencies more precisely than 



E. Checking for systematic errors 

The real experiment is an approximation to the ideal 
scenario of a small-amplitude oscillation in a perfectly 
harmonic trap, as described in the theoretical approach 
in Sec. II. Here we investigate in how far our data are 
influenced by anharmonicities of the trapping potential 
and nonlinear effects arising from the finite amplitude of 
the mode, and we identify the conditions that ensure a 
reliable comparison between the measurements and the 
theoretical predictions. 

The axial compression mode (k = 1) serves us as 
a benchmark to rule out a significant effect of anhar- 
monicities. This mode has been studied extensively be- 
fore [16, 17] and, in the unitarity limit, its frequency is 
temperature independent as long as the gas remains hy- 
drodynamic [23]. Only for very shallow traps, we observe 
a k = 1 mode frequency that is lower than expected. We 
find that trap depths of twice the non-interacting Fermi 
energy are sufficient to obse rve a frequency very close to 
the ideal value of u)k=i = \/12/5uj z throughout the full 
temperature range explored in the present work. Devia- 
tions from this value remain below 0.3% and no signifi- 
cant temperature dependence is observed. 
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FIG. 2. (Color online) Observed collective modes in the time-domain and in frequency space. The three columns refer to 
the mode orders k — 1,2, and 3. We show the normalized density variation signals Sni(z, t)/ni(0) (top row), their Fourier 
transforms Shi(z,u) (middle row), and the mode profiles Sh 1 (bottom row). The measurements were taken for our coldest 
samples with T/Tf fa 0.1. The arrows in the bottom panel show the positions of the repulsive excitation beam for each mode. 
Note that for efficient excitation we adjust the width of the excitation beam to match the local profile at the chosen antinode. 
While, for k = 1 and 3, the beam is centered and addresses the central antinode, the beam is spatially offset for the k = 2 
mode. 



We checked for a possible nonlinear behavior by delib- 
erately overdriving the collective modes. We measured 
the frequency and damping time of each collective mode 
versus the power of the excitation beam. Figure 4 shows 
the results of such a measurement for the k = 3 mode for 
the coldest samples used in our work. The measurements 
show that the frequency stays constant within the mea- 
surement uncertainties up to a power P$ m 1 mW. The 
fact that the spatial profile is already strongly affected 
(see insets) shows that the mode frequency is rather ro- 
bust against nonlinearitics. The mode damping time ex- 
hibits a similarly robust behavior with a slight trend to 
be affected already at somewhat smaller excitation am- 
plitudes. The excitation power Pq is always kept small 
enough to avoid significant effects on the mode frequen- 
cies. 



IV. EXPERIMENTAL RESULTS 

Here we present our main experimental results and 
compare them with the predictions of Sec. II. We consider 
the two modes with k = 2 and k = 3 and discuss how 
their eigenfrequencies and the corresponding mode pro- 
files depend on the temperature. We furthermore present 
data on the temperature-dependent damping of the two 
modes. 

In Fig. 5(a) and (b), we show the measured mode fre- 
quencies Wfe=2 and u>k=3 versus temperature, normalized 
to the axial trap frequency u z . The two limiting cases 
of a T = supcrfluid and a classical collisionally hydro- 
dynamic gas are indicated by the upper and the lower 
horizontal dashed lines in both panels, see Eqs. (2) and 
(3). The theoretical predictions according to Eqs. (4) 
and (5) and the EOS from Ref. [10] are shown by the 
solid lines. For comparison, the hypothetical frequencies 
calculated with the EOS of a non-interacting Fermi gas 
are shown by the dashed lines. For the k = 2 mode we 
have applied both thermometry methods as described in 
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FIG. 3. Evolution of the k = 1 (top), k = 2 (middle) and 
k = 3 (bottom) mode amplitudes at T/Tf ~ 0.1. The black 
circles represent the experimental data for A(t) normalized 
to the maximum value A ma x of |^4(t)|. The red solid lines 
are the fits to the data based on simple damped harmonic 
oscillations. 
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FIG. 4. The normalized oscillation frequency (a) and 1/e 
damping time (b) of the k — 3 mode versus the power Po of 
the green excitation beam [34]. The three insets in (a) show 
the axial density profiles of the atom cloud at the turning 
points of the mode oscillation for three different excitation 
powers. The error bars denote the standard errors obtained 
from fitting A(t). 



the Appendix, with the open symbols representing the 
results from the wing-fit method (Appendix Sec. 2) and 
the filled symbols from the potential-energy method (Ap- 
pendix Sec. 3). For the k = 3 mode we have applied only 
the first method. 

For the k = 2 mode (Fig. 5(a)), the measured mode 
frequencies are in almost perfect agreement with the the- 
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FIG. 5. (Color online) Comparison between experimental and 
theoretical first-sound frequencies for (a) the k — 2 and (b) 
the k = 3 mode. In (a), the reduced temperature T/Tp 
is obtained by the two different methods described in the 
Appendix in Sees. 2 (open red symbols) and 3 (filled black 
symbols). In (b), T/Tf is extracted by only using the first 
method. The theoretical curves (solid lines) are obtained with 
Eqs. (4) and (5) using the EOS of Ref. [10]. For comparison, 
we also show the mode frequencies (dashed curves) that would 
result from the same equations but using the EOS of the ideal 
Fermi gas. In both panels, the upper and lower thin horizontal 
dashed lines mark the zero-T superfluid limits and the clas- 
sical hydrodynamic limits whose values are given by Eqs. (2) 
and (3), respectively. The dash-dot vertical lines in (a) and 
(b) indicate the critical temperature T c /Tf — 0.223(15). 



oretical predictions based on the EOS from [10]. In 
comparison, the disagreement with the dependence that 
would result from the EOS of the ideal Fermi gas high- 
lights the important role of interactions. At the lowest 
temperature realized in our experiment (T/Tp ~ 0.1), 
the frequency lies close to the T = superfluid limit 
(wfc=2/Wz = 2.049), but already shows a significant down- 
shift amounting to almost 1%. This illustrates the high 
sensitivity of the mode frequency to finite-temperature 
effects. At the highest temperatures (T/Tf ~ 0.45) our 
data show a clear trend to go below the asymptotic high- 
temperature value (w/c = 2/^2 = 1.949), which corresponds 
to the classical hydrodynamic case. This nonmonotonic 
temperature dependence can be understood based on the 
first-order correction to the EOS resulting from the virial 
expansion [37] at high temperatures. 

For the k — 3 mode (Fig. 5(b)), the general behav- 
ior is very similar to the k = 2 mode, with the main 
difference that the relative frequency change from super- 
fluid to collisional hydrodynamics [ui]~ = 2/uj z — 2.530 and 
2.280, respectively) is about two times larger. For tern- 



peratures below 0.2 Tp we find similarly good agreement 
as in the k = 2 case. However, for higher temperatures 
there is a significant trend to lie above the predicted fre- 
quencies, so that we never observe values below the high- 
temperature limit. This discrepancy is most likely due 
to dissipative effects, which manifest themselves in mode 
damping. Landau's equations used as a basis for describ- 
ing the collective modes do not contain dissipative terms 
and can therefore not make any predictions on the damp- 
ing behavior. We can therefore address the question of 
damping only from the experimental side. 

In Fig. 6, we show the damping rates of the k = 2 (blue 
open circles) and the k = 3 (red filled diamonds) modes, 
measured at various temperatures. The damping rate 
of the k = 3 mode is always several times higher than 
that of k = 2, and it strongly increases as the cloud gets 
hotter. It is evident that the situation, where we have 
observed significant deviations in the mode frequencies 
(k = 3 and T /Tp > 0.2), coincides with the case of high- 
est damping rates (exceeding 10s _1 ). It is known that 
violation of the ID condition assumed in Sec. II would re- 
sult in a damping of the collective oscillation. Therefore, 
a higher damping rate could naturally be accompanied 
with a larger deviation from the prediction using the ID 
formulation of the Landau's two-fluid model. This sup- 
ports our interpretation of the observed frequency devi- 
ation in terms of dissipative effects. 

For an accurate determination of temperature- 
dependent frequency shifts in our experiment, the k = 2 
mode turns out to be superior to the k — 3 mode. Our 
results indicate that the advantage of the latter mode to 
exhibit larger frequency changes is overcompensated by 
the larger damping, which introduces larger statistical 
uncertainties (see error bars in Fig. 5) and apparently 
also systematic errors. However, for larger trap aspect 
ratios than applied in our experiment, the situation may 
be different and higher modes may provide further inter- 
esting information. 

In Fig. 7(a) and (b), we finally show the observed spa- 




FIG. 7. (Color online) Comparison of experimental mode pro- 
files (data points) and theoretical predictions (solid lines) for 
two different temperatures. In (a), we show the cloud profiles 
obtained with a 2-ms TOF. The solid lines show the density 
profiles obtained from the EOS [10] with T/T F = 0.10 and 
0.45. In (b), we show the experimental and theoretical k — 2 
mode profiles at the indicated temperatures. The compari- 
son for the k = 3 mode is presented in (c). The z scale is 
normalized to the Thomas-Fermi radius Ztf of the zero-T 
interacting gas, which is a factor of f ' ~ 0.78 [10] smaller 
than in the non-interacting case. 
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FIG. 6. (Color online) Measured damping rates of the k — 2 
(blue open circles) and the k = 3 (red filled diamonds) modes 
versus temperature. 



tial profiles of the k = 2 and k = 3 modes, in comparison 
to the corresponding theoretical predictions. For both 
modes we present data sets for the lowest temperature 
that we could realize (T/Tp ps 0.1) and for the highest 
temperature explored (T/Tp m 0.45). In the first case, 
the situation is deep in the supcrfluid regime, whereas the 
second case corresponds to the classical hydrodynamic 
case. For reference, Fig. 7(a) shows the corresponding 
spatial profiles of the unperturbed cloud, from which we 
obtained the temperature following the method of Sec. 3 
in the Appendix. The agreement between the experi- 
mentally observed mode profiles and the theoretical pre- 
dictions is remarkable. Within the experimental uncer- 
tainties and with the mode amplitude being the only fit 
parameter, we find a perfect match. This again highlights 
the validity of the theoretical ID framework and power 
of our experimental approach to higher-nodal collective 
modes. 



CONCLUSIONS AND OUTLOOK 



Appendix: Temperature determination 



We have presented an efficient tool-box to excite and 
detect higher-nodal axial collective modes in a resonantly 
interacting Fermi gas. Our results (see also [24]) reveal 
the pronounced temperature dependence below and near 
the superfiuid phase transition, which is theoretically pre- 
dicted in the framework of a ID two-fluid hydrodynamic 
model [25]. The observed temperature dependence is a 
unique feature of higher-nodal modes and has not been 
observed in any other collective mode studied in Fermi 
gases so far. The excellent agreement of the experimen- 
tally observed mode frequencies with the theoretical pre- 
dictions provides a stringent test for the validity of this 
ID approach and provides an independent confirmation 
of the recently measured EOS [10] of the resonantly in- 
teracting Fermi gas. 

We have also reported first studies on the mode damp- 
ing behavior, which show a strong increase of the mea- 
sured damping rates with the order of the mode investi- 
gated. Dedicated experiments on damping could provide 
valuable information on the viscosity and the thermal 
conductivity of the strongly interacting Fermi gas, which 
may provide further insight into fundamental questions 
related to viscosity [38, 39]. A better understanding of 
damping would also be important to understand the lim- 
itations of the theoretical approach [25] applied to de- 
scribe the modes. 

Generalizations of our experiments to first-sound col- 
lective modes in the BEC-BCS crossover regime [4-6]) 
and to spin-polarized Fermi gases [17, 40, 41] will be 
rather straightforward. A very exciting prospect is the 
extension to second-sound modes [25, 28], where the su- 
perfiuid and the normal component oscillate in oppo- 
site phase. A recent experiment [26] shows the prop- 
agation of second-sound pulses along the trap axis, in 
agreement with a theoretical description based on the 
same approach that is used in the present work. This 
observation may, in principle, be interpreted in terms of 
a superposition of several second-sound modes, but the 
selective excitation and observation of individual modes 
of this kind remains an experimental challenge for future 
experiments. 
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The recorded in-situ density profiles ni(z) arc influ- 
enced by imperfections in the imaging process. While, 
under our experimental conditions, the limited resolu- 
tion and optical aberrations do not play any significant 
role, we have identified another problem (Sec. 1) that can 
considerably affect our thermometry. We here discuss our 
strategies to circumvent this problem, presenting our two 
methods (Sees. 2 and 3) to extract the temperature from 
the observed profiles. 



1. Imperfections of absorption imaging 

In-situ absorption images and images taken with a 
short time of flight (TOF) reveal an apparent reduction 
of the effective absorption cross section, which predom- 
inantly occurs in the denser regions of the cloud. To 
illustrate this effect, we show in Fig. 8 how the apparent 
atom number, i.e. the atom number obtained under the 
assumption of the full absorption cross section, depends 
on the TOF after release from the trap for two experimen- 
tal settings corresponding to a number of N = 1.2 x 10 5 
atoms (red diamonds) and 4.8 x 10 5 atoms (black circles) 
in the trap. Only after a TOF of 2 ms the apparent atom 
number reaches a constant maximum value, which corre- 
sponds to the true atom number. It is evident that our 
in-situ imaging underestimates the actual atom numbers 
by about 15% for the data set with N = 1.2 x 10 5 , and 
by about 30% for N = 4.8 x 10 5 . 

We do not fully understand the reason for this reduc- 
tion, which is clearly related to the high density of the 
cloud, but it is not related to a 'black out' effect, in which 
the imaging light is completely absorbed by a very dense 
cloud. In our case, the maximum optical density is about 
one, which can be properly accounted for by the exponen- 
tial decrease of the transmission with increasing column 
density. We speculate that the underlying mechanism is 
related to multiple scattering of a photon by a few neigh- 
boring atoms when the atoms are very closely packed, 
thereby reducing the effective absorption cross section 
of the atoms. For our geometry, this effect should be 
important when the photon's 'mean free path' becomes 
comparable to the radial size of the cloud, which is indeed 
the case for our experimental conditions. 

By analyzing profiles obtained for different times of 
flight we found that the imaging problem mainly af- 
fects the center of the cloud, where one finds the largest 
density. Therefore the problem does not only result in 
smaller apparent atom numbers, but it also distorts the 
density distribution ni(z) of the cloud as obtained by in- 
situ absorption imaging. To determine the temperature 
of the cloud in spite of this problem, we have adopted the 
two different methods described below. Each method has 
its own systematic uncertainties. Their accuracies can 
only be judged a posteriori by their mutual agreement 
and also by their agreement with certain theoretical pre- 
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FIG. 8. (Color online) Apparent atom numbers obtained with 
different times of flight while keeping other experimental pa- 
rameters unchanged. The red diamonds and the black circles 
represent the data sets with N — 1.2 X 10 5 and 4.8 x 10 5 
atoms, respectively. 



dictions. 

Both methods rely on the accuracy of the total atom 
number TV. We obtain this number through a very care- 
ful calibration of the imaging process and the imaging 
system. For imaging, we typically adopt a 10-^s-long 
imaging pulse and keep the probe intensity less than 4% 
of the saturation intensity. These parameters keep the to- 
tal number of photon scattering events per atom small, 
thereby minimizing the photon-recoil-induced Doppler- 
detuning effect [42, 43] on the light 6 Li atoms. We esti- 
mate the uncertainty in the resulting atom number to be 
typically about 10%, in any case below 20%. 



2. Wing-fitting method 

In the first method, we obtain the cloud's temperature 
by fitting Eq. (9) to the wings of a ID density profile 
taken after a 600 fis TOF [35] , with TV determined inde- 
pendently from a corresponding 2-ms TOF image. The 
wings, where n\(z) is relatively small and the total ab- 
sorption is weak, are essentially free of the image distor- 
tion as described above. An example for a typical wing 
fit is shown in Fig. 9. The Figure also illustrates the dif- 
ference between the true density profile as reconstructed 
by the wing fit (solid line) and the distorted observed one 
(data points). 

In practice, we adjust the number of data points used 
for the wing fit until we obtain an atom number very 
close to the actual value obtained using a sufficiently long 
TOF. Sometimes, we have to use different numbers of 
data points in each wing when the profile on the wings 
is not fully symmetric. This issue is caused mainly by 
optical aberrations in the imaging system. Nevertheless, 
we always find the same T/Tp within an uncertainty of 
10 to 15% using this procedure for samples prepared by 
the same experimental sequence. We also do not observe 
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FIG. 9. Typical fit to the wings of the ID density distribution 
to determine T/Tf from Eq. (9). The filled circles represent 
the data used for fitting, and the solid line is the full profile 
according to the fit. The open circles represent the remaining 
data that were not used in the fit because of the problem to 
determine the correct linear density in the center of the cloud. 
Here a TOF of 600 fis was applied, the total atom number is 
TV = 4.2 x 10 5 , and the temperature is T/T F = 0.36. 



a clear deviation in the so obtained temperature when 
varying the TOF from to 600 /is. Our experience shows 
that the accuracy of T/Tp obtained from this method is 
not very sensitive to slight distortions in the wing profiles. 
Instead, it depends more crucially on the correct atom 
number TV. 



3. Potential-energy method 

In the second method, we characterize the temperature 
in a model-independent way that does not require any a 
priori knowledge of the EOS. Based on the virial theo- 
rem [44] we can obtain the total energy E of the harmon- 
ically trapped cloud from its potential energy, which can 
be calculated from the ID density profile n\(z) according 
to E = 3mu!^ J_ ni(z)z 2 dz. We define a dimensionlcss 
parameter E/Eq, where E is normalized to the energy 
of a noninteracting, zero-temperature Fermi gas with the 
same number of atoms, Eq = |TVfcsT. 

To obtain an accurate value for E/Eq it is essential 
to have accurate knowledge of ni(z). Here we overcome 
the above-discussed distortion problem by 'reconstruct- 
ing' the correct profile. We measure the density pro- 
files of the cloud for two different times of flight, 600-/is 
and 2 ms. The corresponding profiles ?ii(z, 600/xs) and 
n\(z,2ms) provide complementary information. While 
the 600-/is profile provides accurate information on the 
wings, it underestimates the central part of the cloud. In 
contrast, the 2-ms profile suffers in the wings from the 
expansion, but gives accurate information on the cen- 
tral part of the cloud. With an appropriate algorithm 
to combine this information we can reconstruct the full 
profile with reasonable accuracy. The reconstructed pro- 
file shows approximately the correct atom number and it 
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maintains the wing distribution, which has large weights 
ex z 2 in the total energy E. 

We can now convert E/E into T/Tp using the EOS 
from Ref. [10]. Using the universal thermodynamic rela- 
tions of a resonantly interacting Fermi gas, one can show 
that E/Eq is related to the xq = /3/xo at the center of the 
trap by 



_B 4(2^) 1 / 2 Af 2 (x ) 

E o ~ 3(3v / 2^) 1 / 3 M 4/3 (x ) : 

while the xq is related to T/Tp by 



(A.l) 



(27T) 



1/2 



T F [24v^7rAf (a;o)] 1/3 ' 



(A.2) 



Figure 10 shows the conversion between E/Eq and T/Tp 
for a resonantly interacting Fermi gas in comparison with 
the ideal non-interacting Fermi gases. For the resonantly 
interacting Fermi gas at T = 0, E/Eq = -y/f • Here £ 
is the Bertsch parameter, for which Ref. [10] gives the 
value £ = 0.376(4). For the ideal Fermi gas at T = 0, 
E/Eq — 1 by definition. 

We finally note that we applied both methods to var- 
ious data sets to check whether they produce consistent 



results. In general we find satisfying agreement with each 
other, as the example of the data set in Fig. 5(a) shows. 
At very low temperatures, the wing-fit method shows a 
trend to give slightly lower values of T (up to ~10%) as 
compared to the potential-energy method. This indicates 
small systematic uncertainties of our methods. 




FIG. 10. Relation between E/Eq and T/Tf for a 3D harmonic 
trap. The solid and the dashed lines correspond to the results 
obtained using the EOS from Ref. [10] and the EOS of the 
non-interacting Fermi gas, respectively. 
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